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Weak gravitational lensing is one of the key probes of cosmology. Cosmic shear surveys aimed at 
measuring the distribution of matter in the universe are currently being carried out (Pan-STARRS) 
or planned for the coming decade (DES, LSST, EUCLID, WFIRST). Crucial to the success of 
these surveys is the control of systematics. In this work a new method to constrain one such 
family of systematics, known as multiplicative bias, is proposed. This method exploits the cross- 
correlation between weak lensing measurements from galaxy surveys and the ones obtained from 
high resolution CMB experiments. This cross-correlation is shown to have the power to break the 
degeneracy between the normalization of the matter power spectrum and the multiplicative bias of 
cosmic shear and to be able to constrain the latter to a few percent. 



PACS numbers: 98.80.-k, 95.30.Sf, 98.70.Vc 

Introduction. Cosmic shear probes the distribution of 
matter in the universe by measuring the distortions it 
induces in the ellipticities of background galaxies. In 
the past decade, this technique has emerged as one of 
the key probes for cosmology. Thanks to large galaxy 
surveys that are currently under way (Pan-STARRS [1 ) 
or in the works (DES , LSST ^, EUCLID, WFIRST) 
weak lensing promises to tightly constrain the large scale 
properties of the universe and to probe the nature of dark 
energy to unprecedented accuracy [4lfl8]. 

Because systematics are one order of magnitude larger 
than the cosmological signal, their control is critical for 
this experimental program. Several different algorithms 
exist to correct for the instrumental distortions of the 
point spread function and (in the case of ground based 
surveys) for the atmospheric seeing |19l [2Qj . As shown 
by the Shear Testing Programme (STEP) [Ulini, while 
the application of these algorithms corrects for most sys- 
tematics, it can introduce bias in the data. Such biases 
have been catalogued into three generic families: multi- 
plicative and additive biases in shear measurements and 
errors in redshift measurements [19-22\ The focus of 
the present work is multiplicative bias. The lack of pre- 
cise knowledge of multiplicative bias leads to a dramatic 
degradation in the accuracy of the cosmological parame- 
ters thus measured [211 [22] : for example, lack of percent- 
level constraints on multiplicative bias can lead to an 
increase in the errors in the value of flm and erg of 100% 
or more [5T]. Quite generally, multiplicative bias is in- 
sidious because it does not show any scale dependence 
that can be exploited to decouple it from the weak lens- 
ing signal. As such, it is completely degenerate with the 
normalization of the matter power spectrum: if the ob- 
served ellipticities were the only information available it 
would always be possible to trade a non-zero multiplica- 
tive bias for a variation in the total matter in the uni- 
verse responsible for the lensing signal. A first solution 
to this problem was proposed in Vallinotto et al. [23], 
where galaxies' size and luminosity information is used 



to constrain this bias. Such a solution, however, relies on 
somehow quantifying the impact of weak lensing on the 
populations of galaxy sizes and luminosities. 

The present work proposes an alternative way to con- 
strain shear multiplicative bias by cross-correlating it 
with the weak lensing of the Cosmic Microwave Back- 
ground (CMB). Intuitively, these two probes measure es- 
sentially the same observable (albeit extending over dif- 
ferent redshift ranges) but using two completely differ- 
ent techniques. Because the systematics affecting these 
two techniques have completely different natures, cross- 
correlating these two signals is an effective way to con- 
strain their impact. 

The CMB Lensing Field. Experiments aimed at mea- 
suring small scale temperature and polarization fluctu- 
ations of the CMB [21] have recently reached the sensi- 
tivity required to reconstruct the effective deflection field 
arising from the dark matter structures present between 
the last scattering surface and the observer. CMB lens- 
ing has first been detected by WMAP [551 [H] and recent 
measurements from ACT have reported the first detec- 
tion of its power spectrum The measurement of 
the lensing of CMB through quadratic optimal estima- 
tors [28-,-3QJ exploits the statistical properties of the pri- 
mary CMB anisotropics. The dominant systematics for 
this measurement are the ones characterizing high resolu- 
tion CMB experiments: emission from unresolved radio 
and dusty star forming galaxies and thermal and kinetic 
Sunyaev-Zeldovich effects [3T]. As such, these system- 
atics are completely uncorrelated from the ones charac- 
terizing cosmic shear measurements, which arc related 
to the treatment of galaxy images (in particular, atmo- 
spheric seeing and the correction of anisotropics in the 
point spread function) and to the measurement of their 
ellipticities. 

The lensing convergence K{h, Xf) measured for a source 
at comoving distance x f along a line-of-sight (los) in the 
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Survey 


Pixel size 


Galaxies per 


Sky coverage 


O's.N 




(sq. arcmin.) 


nf"" pixel 


(sq. deg.) 




DES 


10 


15 


5000 


0.3 


LSST 


10 


100 


20000 


0.3 



TABLE I: Parameters assumed in the calculation for the two 
cosmic shear surveys. (Ts_jv is the rms of Ks in the absence of 
signal, due to shape noise and measurement errors for a single 
galaxy. 



Survey 


Ic 


Sky coverage 






(sq. deg.) 


Planck 


100 


All sky 


WidePol 


300 


4000 


CMBPol 


1000 


All sky 



TABLE II: Parameters assumed in the calculation for the 
CMB lensing surveys. 
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is defined as 

C r dxWUx,x/-^: 
Jo a{x) 



(1) 



where C = 3n^H^/{2c^) and Wl{x, Xf) = x{Xf-x)/Xf 
is the lensing window function. High resolution CMB 
experiments allow the reconstruction of the convergence 
field K°^^{n, xlss) = i^c^^{fi-)i which extends all the way 
to the last scattering surface. In general, the convergence 
measured from the CMB can be modeled as 



^(n) = 4™°(n)+7i,(n), 



(2) 



where an additive bias nc{n) (or more precisely a set of 
them |35j ) can be introduced by the reconstruction pro- 
cedure |36| . On the other hand, the optimal quadratic 
estimators used to reconstruct the deflection fleld can be 
properly normalized |35j . so as not to introduce multi- 
plicative bias in the CMB lensing signal. In what follows 
a cutoff on the angular modes contributing to the CMB 
lensing signal is introduced. This cutoff is determined so 
that only multipoles that can be measured with signal- 
to-noise ratio greater than unity are included. In this 
(imaging) regime the noise term Uc can be neglected. 

Cosmic Shear and Multiplicative Bias. Weak lensing 
surveys measure cosmic shear in redshift bins of flnite 
thickness. The convergence they aim at measuring is 



\n,Zo) 



dxriix) i^{n,x), 



(3) 



where zq denotes the center of a redshift bin of thickness 
Az and i](x) is a selection function for the given redshift 
bin, normalized so that dxTi{x) = 1- The convergence 
value measured from cosmic shear data in a given redshift 
bin will in general differ from the "true" value because of 
the (possibly redshift dependent) multiplicative bias, so 
that 



KTin,Zo)^b{zo)^r{n,zo). 



(4) 



In what follows a single redshift slice with z g [0.9, 1] is 
considered and on such redshift slice the multiplicative 
bias is assumed to be constant. It is straightforward to 
extend this analysis to other redshift ranges. 



Correlation calculation. Given a pair of surveys, the 
data set will consist of the observed {k°^^,k°^^} over 
the patch of sky where the two surveys overlap. From 
Eqs. (l]|3l, it is straightforward to obtain expressions for 
the elements appearing in the covariance matrix of the 



joint {k 



obs ^obs 



} data set. Let fii and 6ij denote the 



direction of the i-th pixel and the angular separation be- 
tween pixels directed along fii and hj respectively. Fur- 
thermore, let a be the matter power spectrum normal- 
ization, so that P{k) = a^V{k), where 'P(fc) denotes the 
shape of the power spectrum. Then, defining 



gix) = 



dxiWL{x,Xi)v{Xi), 



(5) 



C{x,o,i) ^ ^MW)e~i'/'Ta'r(^^,xy{6) 

and using Limber's approximation, the correlations be- 
tween the different data sets (denoted for brevity by 
{K°^%ni)K°^''{hj)) = {K^K^)ij with {^l,v} = {c, s}) are 
given by the following expressions 



C'b dx 
Jo 



ff(x) W^l(x,Xlss) 



a2(x) 



(7) 
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Cix, 0,j,lc), (8) 



aHx) 



(9) 



where = l^lv/{P^ + ItY^"^- Here Os^n is the rms of 
the shear measurement in absence of signal - due to the 
intrinsic ellipticities of galaxies and measurement errors 
(scaling as iV"^/^, where N is the average number of 
galaxies in a pixel) - and Is is defined as the limiting 
multipole corresponding to the Ks pixel size. Similarly, 
Ic represents the cutoff on the modes contributing to k^. 
As previously mentioned, Ic is fixed to be the limiting 
multipole where the signal-to-noise ratio is equal to one. 

Constraining Multiplicative Bias. The Fisher informa- 
tion matrix is used to estimate the constraining power of 
coupling the Kg and data sets. For the Kg data set, 
estimates are provided for two surveys: DES and LSST. 
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FIG. 1: Projected constraints on multiplicative biases and on the normalization of the matter power spectrum a obtained by 
cross-correlating cosmic shear data from DES (left panel) and LSST (right panel) with CMB lensing data. The red dashed 
curves show results for cross-correlating with Planck data and include a 4% prior on the value of a. The solid black curve shows 
results for DES-|-WidePol (left) or LSST+CMBPol (right) cross-correlations. In the former case, complete overlap between the 
experiments' footprints is assumed. 



For the Kc data set, three different surveys are consid- 
ered: an all-sky survey with low resolution reconstruction 
(Planck |37], Ic = 100), a futuristic all-sky CMB polariza- 
tion experiment leading to high resolution reconstruction 
(CMBPol Zc = 1000) and a ground based polariza- 
tion survey ~ analogous to ACTPol-Wide [39] and SPT- 
Pol [101 - covering 4000 sq. deg. at medium resolution 
(WidePol, Ic = 300). The relevant parameters are sum- 
marized in Tables HI and HIl 

Given a pair of surveys, the data sets will consist of the 
M values {nf^f} and of the N values{K°i;;''}. In general, 
the Fisher matrix is 



(10) 



where A and (p run over the parameters {a, b}, the trace 
is over all the M + N observables and C is the data 
covariance matrix with elements given by Eqs. Be- 
cause of the very large number of k°^^ pixels (1.8 x 10^ 
for DES and 7.2 x 10*^ for LSST), the evaluation of the 
Fisher matrix poses a numerical challenge. It is however 
possible to obtain a reasonable (and conservative) esti- 
mate of the Fisher matrix by noting that all correlation 
functions decrease significantly with increasing angular 
separation. In what follows regions of (3.6)^ sq. deg. are 
treated as statistically independent: over this range of 
separations, the correlations between the different pixels 
drop by more than two orders of magnitude. 

The left panel of Fig. [l] shows the constraints on {a,b} 
projected for DES-fPlanck (red dashed contour) and for 
DES-I- WidePol (black solid contour). In general, cross- 
correlating cosmic shear and CMB lensing allows to break 
the degeneracy between the normalization of the matter 



power spectrum and cosmic shear's multiplicative bias 
and to constrain the latter. Because of the limited num- 
ber of K°^^ pixels that can be reconstructed using Planck 
data (assuming the conservative value of Ic — 100), for 
DES+Planck and LSST+Planck a 4% prior on a (con- 
sistent with current constraints on as) is assumed to 
improve the constraints. In the DES+WidePol case, 
on the other hand, complete overlap between the ex- 
periments' footprints is assumed. Under these assump- 
tions the Fisher matrix estimates show that DES+Planck 
(DES-t- WidePol) data should allow to constrain shear 
multiplicative bias b to about 4% (2%). The left panel of 
Fig.[l]can be directly compared to Fig. 2 of showing 
that this cross-correlation represents a viable alternative 
to using sizes and luminosity information, free from ef- 
fective parameters modeling the impact of lensing on the 
size and luminosity populations. 

The right panel of Fig. [l] shows projections for 
LSST-f Planck (red dashed contour) and LSST-l-CMBPol 
(black solid contour). These data sets benefit from the 
increase in the cosmic shear footprint and the constraints 
on b are respectively reduced to 1.7% and 0.3%. 

Discussion and conclusions. The interplay between 
CMB lensing and cosmic shear data sets allows to lift 
the degeneracy between the power spectrum normaliza- 
tion and cosmic shear's multiplicative bias and to con- 
strain the latter. This is primarily due to the fact that 
the CMB lensing kernel is wide and even if it peaks at 
deeper redshifts, it is broad enough to give a significant 
non-zero correlation with cosmic shear measurements. 

The results reported thus far are based on the quite re- 
markable fact that the CMB lensing signal reconstructed 
from optimal quadratic estimators can be affected by ad- 
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ditive bias but not by a scale independent multiplicative 
bias. In other words, as long as optimal quadratic es- 
timators are used, the CMB lensing signal is properly 
normalized [35] . It is however of practical importance to 
quantify the impact of relaxing such a condition. To do 
this, the CMB lensing multiplicative bias c is introduced, 
so that 



Kf''(7i) = CK*™"(n) + ne(n). 



(11) 



The resulting enlarged set of multiplicative parameters 
{a, b, c} is degenerate with respect to the cosmic shear 
and the CMB lensing data sets. It is nonetheless possible 
to estimate the constraining power of the latter with re- 
spect to b and c by imposing a prior on the normalization 
of the power spectrum. The results obtained by assum- 
ing a 4% prior on a (consistent with current constraints 



on (jg) are reported in Table III Not surprisingly, the 
results obtained strongly depend on the prior assumed. 
However, they also show that under this assumption the 
data allow for testing and constraining both multiplica- 
tive biases at the few percent level. 

It is furthermore possible to speculate that the cross- 
correlation between these data sets could also be ex- 
ploited to constrain the additive biases present in cosmic 
shear and CMB lensing measurements: while CMB lens- 
ing and cosmic shear are characterized by additive biases, 
their cross-correlation is free from these contaminations. 
These developments are left for future investigation. 

Finally, it is possible to remark that the strong syn- 
ergy between these different probes of the dark matter 
distribution suggests that future CMB lensing and cos- 
mic shear surveys should greatly benefit from sharing the 
same footprint on the sky. This fact is not restricted to 
cosmic shear surveys but it applies in general to most 
astrophysical surveys (Lyman-a [4T | |42 ] . BAO, 21-cm): 
because CMB lensing is only sensitive to the dark mat- 
ter distribution, free from any biasing relation, the cross- 
correlation of its signal with other astrophysical observ- 
ables allows to extract biasing relations and to control 
for systematics. 
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